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ABSTRACT 

We test how well available stellar population models can reproduce observed 
^i,^,r,i,z-band photometry of the local galaxy population (0.02 ^ 2; ^ 0.03) as probed 
by the SDSS. Our study is conducted from the perspective of a user of the models, who 
has observational data in hand and seeks to convert them into physical quantities. Stel- 
lar population models for galaxies are created by synthesizing star formations histories 
and chemical enrichments using single stellar populations from several groups (Star- 
burst99, GALAXEV, Maraston2005, GALEV). The role of dust is addressed through 
a simplistic, but observationally motivated, dust model that couples the amplitude of 
the extinction to the star formation history, metallicity and the viewing angle. More- 
over, the influence of emission lines is considered (for the subset of models for which 
this component is included). The performance of the models is investigated by: 1) 
comparing their prediction with the observed galaxy population in the SDSS using 
the {u-g)-{r-i) and {g-r)-{i-z) color planes, 2) comparing predicted stellar mass and 
luminosity weighted ages and metallicities, specific star formation rates, mass to light 
ratios and total extinctions with literature values from studies based on spectroscopy. 
Strong differences between the various models are seen with several models occupying 
regions in the color-color diagrams where no galaxies are observed. We would therefore 
like to emphasize the importance of the choice of model. Using our preferred model 
we find that the star formation history, metallicity and also dust content can be con- 
strained over a large part of the parameter space through the use of ^i,^,r,i,z-band 
photometry. However, strong local degeneracies are present due to overlap of models 
with high and low extinction in certain parts of color space. 

Key words: galaxies: formation - galaxies: evolution - galaxies: stellar content - ISM: 
dust, extinction 



1 INTRODUCTION 

Optical broad band colors have proven to be a powerful tool 
in studying galaxies. Their dependence on luminosity and 
environment have greatly increased our knowledge of these 
systems (Visvanathan Sandage^ 1977 Park et al. 2007 



Lisker et al.|2008 ). Colors are also used to derive quantities 
such as star formation histories and stellar masses (Tinsley 
19681 ISearle et al.|[r 973; Char iot Bruzual||1991 | |Bell 
de jQiig||2001| [Bruzual & Chariot 20031 [Blanton Roweis 
2007| , both of which are key properties for understanding 
galaxy formation and evolution. For example, strong correla- 
tions between stellar mass and galaxy structure (jKauffmann 
et al. 2003), star formation history (Panter et al. 2007), 



chemical enrichment (Panter et al. 2008) and gas content 
[2009 ) have been presented suggesting that 



(Zhang et al 



mass is the main driver behind galaxy evolution. These re- 
lations reflect the importance of gravity on galactic scales 
and, moreover, provide further evidence concerning the ex- 
pected connection between stellar mass and dark matter (cf. 
Moster et al.||2010|. The derivation of star formation histo- 



ries and stellar masses are made through a modelling of the 
light emission from the galaxy and, if necessary, modelling 
of obscuration by dust. The method enables a comparison 
between observational quantities and galaxy formation mod- 
els (e.g. De Lucia et al.|[20Q6[ Bower et al. 2006; Guo et al 
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|2011| ). The quality of the derivation of star formation his- 
tories and stellar masses from colors naturally depends on 
the quality of the models. A straightforward test is to check 
how well the models can reproduce the ensemble of observ- 
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ables. In this paper we therefore take a closer look at how 
successful various models are in reproducing i^,^,r,i,z-band 
photometry of the /ocaQ galaxy population. 

The base ingredient of stellar population models of 
galaxies are single stellar populations (SSPs), which are 
combined into star formation histories by linear combina- 
tions (Tinsley 1968 Searle et al. |1973J . SSPs can conse- 



quently be seen as a basic ingredient of the models having 
a great impact on the emergent spectral energy distribu- 
tion. A large number of single stellar population models are 



available in the literature, including those of Leitherer et al, 



( 1999| ) ;|Bruzual k Chariot] ( |2Q03t ; |Marastont (2005 J ; ,Kotulla 



et al. (2009), which we use in this paper. We do not aim at 



being complete considering the numerous options available. 
However, we intend to cover some of the most widely used 
models predicting i^,^,r,i,z-band photometry. 

We treat dust obscuration in a simple way. The am- 
plitude of the extinction is coupled to a galaxy's star for- 
mation history, metallicity and viewing angle, as motivated 



by works of Cid Fernandes et al. |2005); Engelbracht et al. 



(2008); da Cunha et al. (2010); Masters et al. (2010), and 



a power law form of the the wavelength dependence of the 
extinction is assumed fsee, e.g., [Chariot F^[200^ . This 



method can be seen as an easy, but simplistic, approach with 
respect to detailed modelling of dust (e.g Tuffs et al.]|2004' 



Popescu et al.|[20TT ), though it has the nice feature that it 



only requires the axis ratio apart from colors. 

As spectroscopic data carry a lot more information than 
i^,^,r,i,z-band photometry a comparison with these kinds 
of data can serve as an important test. We therefore carry 
out a detailed comparison between literature data based on 
spectroscopy and our photometric model in order to validate 
the latter. 

We note that a different approach to the testing of stel- 
lar population models has been performed by |Conroy et al.] 
(2009[|2010|) andjConroy Gunn| ([20To| . We refer the reader 



to this series of papers if they are interested in more details 
on the impact and uncertainties of various model ingredi- 
ents rather than overall performance for optical broad band 
photometry. 

This paper is organized as follows. In Section |2] the 
observational sample used for the model evaluation is pre- 
sented. The SSPs are introduced in Section |3] and how these 
are combined into star formation histories is described in 
Section [4] How we treat dust is described in Section [5] The 
performance of the models is presented in Section [6] The re- 
sults are discussed in Section [7] and summarized in Section 

El 



2 OBSERVATIONAL SAMPLE 

Our source of observational imaging data is the Sloan Digi- 



tal Sky Survey (SDSS) data release 7 (DR7 Abazajian et al. 
[200 9 ). Galaxies with measured spectroscopic redshifts from 
the SDSS in the range 0.02 < z < 0.03 are chosen to mini- 
mize the shifts in the spectral energy distribution caused by 
redshifts, while keeping a large sample with photometry of 



sufficient quality (note that, in particular, nearby galaxies 
have inaccurate photometric measurements from the SDSS 



pipeline ( Blanton et al. 2005 [Lauer et al. 2007 Bernardi 



et al.||2007| [Lisker et al.||2007| [Blanton et al.||2011t .) Our 

sample consists of 24120 galaxies. We downloaded model 
magnitudes (modelMag) for these galaxies, which have been 
corrected for Galactic extinction using the maps of ISchlegel) 
et al. ( 1998 ). The sample spans an absolute magnitude range 
of —17 ^ Mr ^ —23, and thus only the brightest of the 
dwarf galaxies are included. The photometry is k-corrected 
using the code of jChilingarian et al.| ( [2010 ) , but the small 
amplitude of the corrections 0.1) suggests that the errors 
introduced in this step are negligible. In this paper we only 
make use of magnitudes in the AB system, including the 
SDSS to AB conversion factor^of uab — usdss — 0.04mag 
and zab = zsdss + 0.02mag. A quality assessment of the 
galaxy colors, as described in the Appendix [A] shows that 
the distribution of colors for the sample appears to be con- 
sistent with independent measurements. 

Information about galaxy structure can aid in the in- 
terpretation of the modelling and is therefore of interest. 
Fig. ^ illustrates the properties of our galaxy sample, rely- 
ing on structural parameters from the SDSS pipeline and the 
Galaxy Zo o project ([Lintott et al.'2008). The SDSS pipeline 
( [Stoughton et al.|2002 ) models galaxy light profiles through 
a best-fitting linear combination of a de Vaucouleurs and an 
exponential profile. The de Vaucouleurs fractions, fracDeV, 
is a measure of the fraction of light in the de Vaucouleurs 
profile of the two component fit. We further make use of 
the ratio of the semi major, a, and semi minor axis, 6, at 
the object's 25mag/arcsec^ isophote which we denote as a/b. 
From the Galaxy Zoo project (iLintott et al. 2011) we ob- 



tain Pell which gives the probability that the objects are as- 
signed elliptical morphology through a visual classification. 
We chose the values that were debiased from resolution ef- 



fects (see Bamford et al. (2009)). Moreover, we make use of 
absolute r-band magnitudes, Mr, computed from the SDSS 
modelMags and redshifts using i7o=71km/s. 

To test the performance of the photometric modelling 
we compare the outcome with results based on spectroscopic 
data from the SDSS. As the spectra are obtained using op- 
tical fibers with a diameter of 3" only a fraction of the light 
from each galaxy is sampled. Stellar mass weighted ages and 
metallicities as well as total IZ-band extinctions. Ay, are 



taken from the SEAGal/STARLIGHT project (Cid Fernan- 



des et al. 2005 ). These data have been derived from the SDSS 



fiber spectra using a spectral synthesis technique in which 
dust is treated as a foreground screen with a Cardelli et al.| 
( 1989| ) extinction law for i?v=3.1. Luminosity weighted ages 
and metallicities are compiled from jGallazzi et al. (2005), 
who derived ages also by a spectral synthesis technique, 
while metallicities are instead derived through an analysis of 
absorpti on line indices, see (.Worthey^l994,; Bruzual Char-| 



lot |2003| [Thomas et al.|2003| ). We moreover make use of spe- 



cific star formation rates based on Brinchmann et al. 



(2004) 



with an improved aperture correction (Salim et al.|2007). All 



data are based on SDSS DR7 except those compiled from 



Gallazzi et al. ( 2005|) which are ba sed on SDSS data release 
4 ( Adelman-McCarthy et al.|2006 ). For convenience we only 



^ By focusing on local galaxies (0.02 < z < 0.03) we circumvent 
the shift in the spectral energy distribution caused by redshift. 



^ See www.sdss.org/dr7/algorithms/fluxcal 
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Figure 1. Properties of the galaxy samples used in this work. 
Top panels: 2D histograms in the (i-z) vs. (g-r) plane of the 
number of galaxies in the full sample selected from SDSS imag- 
ing data (left) and the spectroscopic subsample (right). Middle 
and bottom panels: Average values of absolute r-band magni- 
tudes, Mr, r-band de Vaucouleur fractions, fracDeV, debiased 
probabilities of having elliptical morphology, Peii, and axis ratios, 
a/6, for the full sample in bins of (i-z) vs. (g-r). Note that Peii is 
only available for about 90% of the full sample. 



use spectroscopic data for galaxies in common between these 
studies. The spectroscopic sample used throughout is there- 
fore smaller than the photometric sample and contains 7542 
galaxies in the range 0.02 < z < 0.03. To allow a fair com- 
parison with properties derived from photometry we make 
use of magnitudes measured within the spectroscopic fibers 
when comparing with the data from |Cid Fernandes et al.| 
( |200 5) and Gallazzi et al."('2005), while the modelMags are 
used for the specific star formation rates, since aperture cor- 
rections based on integrated SDSS photometry have been 
applied to these. 



3 SINGLE STELLAR POPULATION MODELS 

The stellar population models we test in this paper are sum- 
marized in Table [T] We retrieve SSPs for each of these mod- 
els which we combine into star formation histories as de- 
scribed in the following section. As most of the models have 
several options regarding specific model ingredients we de- 
cided to constrain the choice to avoid ending up with an 
impractical number of models. Details of the model sub- 
selection are given in Table [l] and the model ingredients 
are briefly commented below. The total amount of spectral 
information contained in i^,^,r,i,2;-band photometry can be 
fully represented in four dimensions by, e.g., the (u-r), (g- 
r), (i-r) and (z-r) colors. However, for practical reasons we 
make use of 2D projections in the form of color-color dia- 



grams, and for simplicity we constrain ourselves to the use of 
the (u-g)-(r-i) and (g-r)-(i-z) planes. Note that these projec- 
tions alone can not capture all available information, though 
they ought to include the majority since they contain infor- 
mation from all five bands. 



3.1 Stellar spectra 

The models fall in two different categories regarding how the 
stellar spectra are generated: theoretical models based on 
stellar atmosphere models ( Lejeune et al.|1997 ) and models 
based on libraries of observed stellar spectra (Le Borgne 
et al.||2003| ). 



3.2 Stellar evolution 

In terms of stellar evolution several different sets of mod- 



els are used, including Schaller et al. (1992); Bertelli et al. 
([1994)- Cassisi et al. (2000);Pirardi et al.| ( |2000J ; ,Marigofc 



Girardi|(2007); Marigo et al. (2008) 



3.3 Initial mass function 



In a recent review Basti an et al.| (|2010') conclude that there is 
no clear evidence that the initial mass function (IMF) varies 
strongly and systematically, and that the majority of sys- 
tems on galactic scales are consistent with having a |Kroupa| 
( 2001| ) or Chabrier ( |2003| initial mass function. However, an 
IMF that evolves with time cannot be excluded and recent 
studies do suggest a steeper IMF in massive elliptical galax- 



ies ( |van Dokkum k Conroy|2010|[2QlT||Spiniello et al.|2011 ). 

Here we investigate the the behavior of our models with re- 
spect to single slope IMFs. The |Leitherer et aL ( |1999| ) model 
(S99) has the option of freely choosing the IMF. Figure |2] 
shows the behavior of the model boundaries as a function 
of IMF slope. Clearly, the S99 models do not reproduce the 
observed distribution of colors very well regardless of IMF 
slope. As expected, the colors only change mildly with the 
IMF and we therefore chose to apply either a [Kroupa| ( |2001| 



Chabrier] ( |2003 ) IMF. For each of the models one can 



choose between one of these two IMF presciptions. 



3.4 Emission lines 

Emission lines are prominent features in optical spectra of 
star forming galaxies. It is therefore desirable to include 
emission lines in the modelling of photometric properties. 
Two of the models we use incorporate gas emission, namely 
S99 and GALEV ( Kotulla et aH2009) . The former predicts 
the strength of Ha and H^^ while the latter predicts several 
other lines in the optical, which are expected to be important 



in star forming galaxies (Anders & Fritze-v. Alvensleben 
2003). Modelling of line emission from SSPs inevitably suf- 



fers from some uncertainty, but can be sufficiently accurate 
for predicting broad-band colors ( Gyory et al.||20lT ). From 
the color evolution of the SSPs it is clear that gas emission 
in the models only has a significant effect on the colors at 
ages of about lO^yr or younger. Fig. [s] shows how the out- 
lines of the model grids shift with the inclusion of emission 
line modelling. Ha and H/3, the only emission lines included 
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Figure 2. Top panel: {i-z) versus {g-r) for our sample of SDSS 
galaxies (number density map) along with model boundary con- 
tours (colored lines) for the S99P models with various IMF slopes 
from 2.1-2.9 (Salpeter=2.35). Bottom panel: Same as above but 
for (r-i) versus {u-g). 



in S99, only affect the r and g bands while additional promi- 
nent emission lines included in GALEV also affect i and 
z causing large discrepancies between the two models in (r- 
i). Modelling of line emission as in S99 or GALEV can be 
included in any of the models we consider by adding the 
corresponding change in the colors of S99 or GALEV. This 
can only be done for SSPs of the same age and metallicity 
under the assumption that the stellar emission in u,g,r,i and 
z is similar in the two models. 



4 STAR FORMATION HISTORY 

We synthesize star formation histories from the SSPs 
through the use of smooth models to which some stochas- 
tic sampling has been added. Such Monte Carlo libraries 
of star formation histories have previously been employed 



by several authors (Kauffmann et al.[2003 da Cunha et al. 



2008| [Zibetti et al.|2009t and have proven versatile in mod- 
elling stellar population properties. The stellar population 
synthesis method we employ is thus not new, but the details 
differ between our and these previous works. The addition 



Figure 3. Top panel: {i-z) versus {g-r) for our sample of SDSS 
galaxies (number density map) along with model boundary con- 
tours for the GALEV (red and yellow lines) and the SB99 (blue 
and green lines) models, with and without the inclusion of emis- 
sion line fluxes. Bottom panel: Same as above but for {r-i) versus 
{u-g). 



of some random component in the modelling is motivated 
by the knowledge that star formation to some extent occurs 



stochastically. Processes such as galaxy merging ( Di Matteo 
et al.|2007| and ram pressure stripping fGunn Gott|1972 ) 
often occur on short timescales and can strongly influence a 
galaxy's star formation history. 



Gavazzi et al. (2002|), inspired by the work of Sandage 
(11986), showed that a "delayed exponential" star formation 



history does a better job in reproducing colors of Virgo clus- 
ter galaxies than the classical exponential model. We there- 
fore use this model as a starting point for our star formation 
histories 



SFR: 



rj-\2 



(1) 



where SFR is the star formation rate, T is the time from the 
initial onset of star formation and ip is a, parameter governing 
the decline of star formation over time. We take the galaxy 
formation time to be 13.5Gyr ago. The exact starting point 
is not crucial due to the slow evolution of spectral properties 
at old ages, but it is to some extent motivated by current 
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Model ID 



reference 



stellar spectra 



stellar evolution 



IMF 



line emission 



BC03hr 

BC031r 

CBOThr 

CBOTlr 

S99P-^ 

S99G-^ 

GALEV^ 

M05 



Bruzual & Chariot (20031 
^Bruzual Chariot^ ( ^2003 1 



m prep. 

in p rep.*^ 

Vazquez & Leitherer ( 2005| > 
LeithereT et al."('l999| 
Kotulla et al. (2009) 
|Maraston| |2005| > 



Le Borgne et a l. J 2003 1 



Lejeu ne etaD(' 19971 
Le^orgne et al. (2003|) 

^eune et al. ( 1997 
Lejeune et al. (1997) 
Lejeune et al. (1997) 
Lejeune et al. (1997) 
[Lejeune et ar| ( [T997 ) 



Bertelli et al. 



Bertelli et al. 



Marigo et a l. 



Marigo et al. 
Uirardi et al. 
Hchaller et al. 



Bertelli et al. 



. et al. 



( 1994 ) 

mm 

(2008)" 
(2008, 
(201 
(1992) 
(1994)' 
|( |2000p 




fchabrier (2003) 
tehabrier (2003) 
Chabrier (2003) 
Chabrier (2003) 
roupa (2001) 
roupa (2001) 
Kroupa (2001) 
|Kroupli| ( pOOTf 



no 
no 
no 
no 

yes 
yes 
yes 
no 



^|fe|c|rf|e|i 



^ Additional references left out, see model reference. 

^ These models have the option of choosing between populations with or without blue horizontal branch. We chose the latter. 
^ S. Chariot private communication. 
d 



Updated SDSS filter responses from 



Doi et al, 



(2010) used. 

GALEV uses the 1999 version of isocrones trom tlie Padova group. 
f Magnitudes obtained by convolving the low resolution spectra with the updated SDSS filter response from 

Table 1. Sources of SSP models 



Doi et al. 



(2010) 



estimates of the onset of reionization in the Universe (Fan 
[eTal 2006). 

Star formation histories are created through a sampling 
of Eq. 1 using SSPs for values of ip in the range IGyr to 
20Gyr. The sampling is done logarithmically, since the color 
evolution is roughly proportional to the logarithm of the age. 
Thus, the age of each SSP, t^, is randomly drawn from the 
following distribution 



exp 



2^2 



(2) 



To compensate for the logarithmic sampling the strength, 
i.e. mass, of each SSP is multiplied by its age. In total 5ip 
SSPs are used for each star formation history (more SSPs 
are needed to sample the larger range in log (age) spanned by 
models with higher values of ip.). Stochasticity is introduced 
in the model through the sampling and by further modify- 
ing the strength of each SSP by multiplying with a random 
component. The random component is drawn from a distri- 
bution which is taken to be the absolute value of a normal 
distribution with (/i = 0,cr = 1). To include many different 
metallicity configurations the SSPs constituting each SFH 
are randomly divided into five groups and each of these 
groups is assigned a random metallicity from the options 
available for that particular SSP. 

The chosen input SSPs limit the lowest age s that can 
be includec j^ However, note that in the model by Chariot 



Fall 



( 2000| the light from populations younger than 10 yr 
is heavily obscured by the dust in the birth clouds and does 
not contribute significantly to the integrated light at optical 
wavelengths. We do therefore not need to consider these 
missing models at young ages. 

Using optical spectral energy distributions the resolu- 
tion in terms of the age of a stellar population seems to be 
limited to at least ^ 10% (Gonzalez Delgado & Cid Fernan- 
des|2010| ). Any gaps introduced by a discrete sampling of a 



^ Note than |Maraston| ( [2005| do not offer young ages at the lowest 
metallicity. 



star formation history can thus safely be neglected if they 
are smaller than about 10%. 

For a range of valueQwe create 1250 realizations of the 
different SFHs, which results in 50000 models for each set 
of SSPs. Contours of the model boundaries in color-color 
space are created in the following way. 2D-histograms in 
{g-r)-{i-z) and (u-g)-(r-i) are made with bin sizes of 0.0055, 
0.0040, 0.0100 and 0.0035 for (g-r), (i-z), (u-g), and (r-i), re- 
spectively. These histograms are convolved with a Gaussian 
kernel (a = 3pixel) to remove the noise at the edges after 
which contours are drawn at a density of 0.2 models/pixej^ 
Figure |4] and |5] shows the outlines of these stellar population 
models (as summarized in Table [T]) in the [u-g)-(r-i) and 
{g-r)-{i-z) planes. 



5 DUST 

In this paper we employ a simplistic model for dust extinc- 
tion in galaxies that has, by construction, no free param- 
eters. Although several assumptions are made we expect 
our model to work reasonably well since it is constructed 
to capture the most important features of dust extinction in 
galaxies. It is motivated by: 

1) The connection between the presence of dust and 



star formation/young stellar populations, da Cunha et al. 
(2010j) showed that the dust mass in galaxies can be esti- 



mated remarkably well using the average star formation rate 
during the last lO^yr (both quantities are determined from 
a model fit to photometry over a large wavelength range, 
^ 0.1 — 100/im). They use a relation in the form of a power 
law with index 1.1 and find a scatter of about 0.5dex over 
at least three orders of magnitude in star formation rate. 
If this relation is normalized by the stellar mass, then the 
mass fraction of the dust is proportional to a power of the 



^ = 1.0, 1.5,. ..3. 0,3. 25, ...,9.0,10.0, ...,20.0 in V^. A variation 
in step size turned out to be useful for better coverage in the 
(u-g)-(r-i) and (g-r)-(i-z) plane s. 

^ The model boundaries in Fig. 2|6 are created in the same way. 
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is essentially a smoother version of the lO^yr cut adopted 




1.2 1.4 1.6 1.( 

(u-g) 



Figure 4. (r-i) versus (u-g) for our sample of SDSS galaxies 
(number density map) along with model boundary contours for 
eight different sets of stellar population models (color coded con- 
tours). Details of the models are summarized in Tableland de- 
scribed in Sect. [3] and [4] 




0.4 0.6 

(g-r) 



Figure 5. {i-z) versus {g-r) for our sample of SDSS galaxies 
(number density map) along with model boundary contours for 
eight different sets of stellar population models (color coded con- 
tours). Details of the models are summarized in Tableland de- 
scribed in Sect. [3] and 2] 



recent specific star formation rate, a property that can be 
constrained by optical colors alone ( [Brinchmann et al.|2004 ). 
We take the effective r-band extinction, Ar, to be 



Ar Aq 



i=0 



-(ti/to) 



i=0 



V 



Jo 



IS.SGyr 



-(*/to)di 



(3) 



/ 



where Aq and to are constants, t^ is the age of the ith SSP, si 
is the strength of the ith SSP, the summation is done over all 
n SSPs in the star formation history and the power index 



of 1.1 is taken directly from da Cunha et al. (2010). This 



equation makes Ar dependent on the fraction of young stars 
by comparing the strengths of the SSPs weighted by the 



factor e 



-(ti/to) 



(numerator) with the corresponding value 



by da Cunha et al. (2010) assuming that Ar is proportional 
to the dust to stellar mass ratio. Aq is constrained by the 
reddest galaxies observed. A value of Aq = 0.40 keeps the 
models within or close to the observed cloud of galaxies, to 
should be chosen in accordance with the work of Ida Cunhal 
et al.|(,2010) and we adopt to=3 • lO^yr. 

2) An observed correlation exists between metallicity 
and dust to gas mass ratio. We base the metallicity depen- 
dence of our dust models on the observational results of 
[Engelbrac ht et al.| (|2008| . Their relation between the nebu- 
lar metallicity, Z, and HI to dust mass fraction (Mni/Mdust) 
can be rather well fitted by a power law (the fit has a scatter 
of 0.5dex of which at least half can be explained by mea- 
surement errors) with low metallicity galaxies having higher 
Mni/Mdust- Assuming that the nebular metallicity is the 
same as the average metallicity of the stellar populations 
contributing to the dust (see Eq.|3|, z, this leads to a mod- 
ification of Ar according to 



Ariz/zeY 



(4) 



where m is a constant. Here we have further assumed that 
the metallicity is the only parameter governing Mni/Mdust 
for a fixed star formation history and that Mni/Mdust is 
proportional to Ar. We use m=1.7 as measured from the 
data presented in Engelbracht et al. ( 2008| ). We note that 



De Lucia & Blaizot (2007) also use a power law to model 



the metallicity dependence of the effective extinction with a 
power law index very similar to the one we adopt. 

3) At optical wavelengths, the wavelength dependence 
of the extinction appears to be rather well modelled by a 
power law of index k ( [Chariot k Fall||2000 ) 

(5) 



Ax oc A'^ 



where Ax is the effective extinction at a wavelength A. | Char-] 
lot & Fall ( 2000 ) use a different proportionality constant for 



stellar populations older and younger than lO^yr. However, 
the contribution of optical light from young populations 
(< lO^yr) is small due to obscuration from the dust clouds 
in which the stars were born. We therefore apply a single 
power law to model the wavelength dependence of the 
extinction. This should be seen as a simple approach, as 
we do not know how the detailed effective extinction curve 
varies over a range of galactic environments. The colors of 
dusty models depend on the constant k. We use A:=-l.l 
which is an appropriate choice considering the width of 
the model cloud in (r-i) at (u-g)>1.3 as compared to the 
observations. The effective extinction in u,g,r,i and z is 
computed using the effective wavelength of these filters. 
The use of effective wavelengths rather than a flux weighted 
average only introduces errors in A of a few percent (as 
estimated using the S99 models without emission lines). 

4) The extinction is dependent on the angle under which 
the galaxy is seen. Empirical relations between colors and 
inclination for spiral galaxies in the SDSS were studied by 
Masters et al. (2010). Here we adopt a similar functional 



for a constant star formation history (denominator). This 



form to express how the effective r-band extinction is mod- 
ified depending on the isophotal axis ratio, a/6, 

= Ar + 7rlogio(a/6) (6) 
where 7r is a constant. To determine 7r we select spi- 
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ral galaxies {Peii < 0.2) from our spectroscopic sam- 
ple and plot Av from |Cid Fernandes eT al. ( 2005 ) ver- 
sus \ogiQ{a/b). For four bins in fracDeVr we determine 
jv (in analogy to jr) through least square fits. The bins 
were chosen such that they contain equal numbers of galax- 
ies and we find {fracDeVr,jv) = (0.44-1.00,0. 51), (0.17- 
0.44,0.58), (0.02-0.17,0.57) and (0.00-0.02,0.58). Given the 
small bin to bin variations we adopt 7r=0.47 (i.e. 7^=0.55) 
and apply Eq. l6] for all galaxies having Ar > 0.05mag ac- 
cording to Eq.[3|and[4] The extinction for galaxies with little 
or no dust are thus assumed to be inclination independent 
while dusty galaxies are treated as spirals. 

We have constructed our model so that it does not have 
any free parameters: The extinction is completely coupled 
to the galaxy's star formation history, chemical enrichment 
and the observed axis ratio. The constants included in the 
model have values that can be inferred from previous stud- 
ies, but some were adjusted in order to produce as realistic 
an ensemble of models as possible as judged by comparing 
the distribution of observed and model colors. We adopted 
Ao 0.40, to O.SOGyr, m 1.7 and k -1.1. 

Note that models which do not have any recent (t < 
to) star formation are essentially unaffected by extinction 
and devoid of line emission. These models are simply linear 
combinations of single stellar populations. 



6 PERFORMANCE 

We test the models described above using an observational 
sample of galaxies from the Sloan Digital Sky Survey DR7 
( [Abazajian et al. 2009') as described in Sect. [2] Given the 
boundaries of the model cloud for dust-free models in Fig. 
|4]and[5] in comparison with the observations we see large 
differences between the various models. In particular the 
models of Chariot and Bruzual in preparation (CB07) and 
GALEV occupy a large region in (g-r)-{i-z) where no galax- 
ies are found to reside and the same is true for CB07 in {u- 
g)-{r-i). The overall best agreement between observations 



and models is shown by the Bruzual & Chariot ( 2003 ) mod- 



els based on observed stellar spectra (BCOShr) and spec- 
tra from stellar atmosphere models (BC031r). Of these two 
BCOShr performs significantly better in the reddest region of 
both color-color diagrams. In the following we have therefore 
adopted BCOShr as our preferred model. A missing ingredi- 
ent in this model is the line emission. As we showed in Fig. 
[S] inclusion of line emission from GALEV only changes the 
model track slightly. We can therefore safely include line 
emission from GALEV into the BCOShr models as discussed 
in Sect. |S.4| GALEV is preferred over SB99 since the for- 
mer takes more emission lines into account. Moreover, we 
apply our dust models to BCOShr. The resulting change in 
the model boundaries is shown in Fig. [6] Considering the 
photometric errors the model performs very well in repro- 
ducing the colors of the galaxy sample, except at about (g- 
r,i-2;) = (0.4,0.0), where models appear to be missing. The 
colors of the SSPs in Fig. [6] show that galaxies dominated 
by recent star formation can have (i-z) <0.0. However, these 
objects would be very blue in (g-r) and thus cannot explain 
the "missing models" . Considering Fig. [6] our recipe for the 
star formation histories combines the blue points into mod- 
els within the green contours, which thereafter are shifted 




(g-r) 




- BC03hr+em+dust 
BCOShr 

SSPs (BC03hr) 



Figure 6. (i-z) versus (g-r), upper panel, and (r-i) versus (u- 
g), lower panel, for our sample of SDSS galaxies (number density 
map) along with model boundary contours for the BCOShr model 
(green) and the BCOShr model with dust and line emission (red). 
The SSPs from BCOShr with ages older than lO^yr are shown 
(blue dots, metallicity increase from light to darker blue) and 
these are connected along lines of constant metallicity. 



to the red contours by dust and line emission. Fig. [5] clearly 
shows that a mixture of SSPs of different ages and metal- 
licit ies is needed to explain the locus of observed galaxies. 
Additionally, dust is needed to explain some of the reddest 
galaxies. 



6.1 Derived parameters 

To illustrate the predictions of the models in terms of galaxy 
properties Fig. tJIq show the stellar mass weightecj^ ages, 
metallicities and total r-band extinction in the color-color 
planes along with their standard deviations. This is shown 
for the BCOShr models with dust and line emission which 
are in good agreement with the observations (Fig. |4] and 
[s]). A strong, rather orthogonal dependence is seen between 



° In this work stellar masses are derived using a |Chabrier| Impost 
IMF and the stellar mass are taken to be the mass in present day 
stars and remnants. 
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Figure 7. Left panels: stellar mass weighted age (color coded 
pixels) as a function of (i-z) versus {g-r)^ as well as {r-i) versus 
{u-g) for our 50000 models based on BCOShr with dust. Right 
panels: standard deviation in stellar mass weighted age (color 
coded pixels) as a function of {i-z) versus {g-r), as well as {r-i) 
versus {u-g) for our 50000 models based on BC03hr with dust. 
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Figure 8. Left panels: stellar mass weighted metallicity (color 
coded pixels) as a function of {i-z) versus {g-r), as well as {r-i) 
versus {u-g) for our 50000 models based on BC03hr with dust. 
Right panels: standard deviation in stellar mass weighted metal- 
licity (color coded pixels) as a function of {i-z) versus {g-r), as 
well as {r-i) versus {u-g) for our 50000 models based on BC03hr 
with dust. 



age and metallicity. However, for a region of color space 
roughly from {u-g,r-i) = {l.l,0.2) to {u-g,r-i)={1.6,0A) and 
{g-r , i-z) = {0.5, 0.2) to {g-r, i-z) = {0.8, 0.3) all properties are 
rather poorly constrained due to a mix of models with high 
and low dust extinction. The model library was constructed 
to cover the {u-g, r-i) and {u-g, r-i) diagrams and does not 
tell anything about the probability that a specific model is 
a good representation of the observations. Further improve- 
ment can be achieved using a Bayesian approach. With help 
from the semi-analytic galaxy formation models of |Guo et ah] 
(2011 ) we compute model weights to match the distribution 



of stellar population models with semi-analytic galaxies in 
terms of density in stellar mass weighted age versus stellar 
mass weighted metallicity, see AppendixjC] For 10 randomly 
chosen example galaxies we show the SDSS color compos- 
ite imag^ along with with a star formation history derived 
through a Bayesian maximum likelihood for the BCOShr 
models with emission lines and dust, see Fig. |10[ 



6.2 Visual checks 

As a test we selected random samples of galaxies for three 
different parts in {g-r)-{i-z) space, which according to our 
models have certain special characteristics. The first group 
contains galaxies that belong to a well populated region 
around {g-r , i-z) = {0.4, 0.0) that can not be reached by any 
models (Fig. 



11 



yellow dots) . This group seems to be popu- 
lated by star forming galaxies that in many cases exhibit Sm 
type morphology. The second group at {g-r , i-z) ^{0.8, 0.3) 
contains galaxies that, according to our model, are expected 
to have strong extinction by dust (Fig. 11 green dots). This 



group contains many edge on or highly inclined disk galaxies. 
The last group at (^-r, 2-^)^(0.8,0.2) is expected to contain 
dust-free galaxies (Fig. |ll| grey dots). With a few exceptions 



Images are taken from http:/ /cas. sdss.org/dr5/en/tools/chart /list 
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Figure 9. Left panels: total r-band extinction (color coded pix- 
els) as a function of {i-z) versus {g-r), as well as {r-i) versus {u-g) 
for our 50000 models based on BCOShr with dust. Right panels: 
standard deviation in total r-band extinction (color coded pixels) 
as a function of {i-z) versus {g-r), as well as {r-i) versus {u-g) for 
our 50000 models based on BCOShr with dust. 



this group contains early-type galaxies with little or no visi- 
ble dust and star formation. SDSS multicolor composites of 
individual representative galaxies are shown in Fig. |11[ 



6.3 Comparison with spectroscopic studies 

There is no question about the superior information content 
of optical spectra with respect to broad band photometry. 
Photometry is nevertheless an important probe of stellar 
populations of galaxies as it is easy to obtain and because 
spectroscopic data often suffer from aperture bias. As men- 
tioned in Sect. [2] we compiled spectroscopically determined 
galaxy parameters from the literature for testing the per- 
.aspformance of our photometric model. The parameters con- 
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Figure 10. Images (left): 50'x50' SDSS u,g,r,i,z color composite images ( [Szalay et al.| 2002) of 10 randomly selected galaxies from our 
sample. Diagrams (middle): (i-z) vs. (g-r) colors for the galaxies (red dots) and the model boundaries (black lines). Diagrams (right): 
Normalized star formation rates (red lines) with errors (shaded region) for the corresponding galaxies as a function of stellar age. These 
star formation histories have been derived through a Bayesian maximum likelihood using all 10 colors and our models based on BCOShr 
with emission lines and dust (see Appendix [c}. Stellar mass weighted ages and metallicities, z, as well as effective r-band optical depths 
are indicated above each image+diagram pair along with absolute r-band magnitudes, Mr, and Ho; equivalent widths from the SDSS 
(positive values indicate emission and negative values absorption). 
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sidered are stellar mass weighted ages, agcmass, luminos- 
ity weighted ages, agcr, stellar mass weighted metallicities, 
Zmass, luminosity weighted metallicities, zium, specific star 
formation rates, 5'5'F/|^ r-band mass to light ratios, M/Lr, 
and y-band extinctions. Ay- To be able to make a fairer 
comparison between photometrically and spectroscopically 



^ Specific star formation rates are taken as the ratio of the av- 
erage star formation rate during the last lO^yr and the galaxy's 
stellar mass. An average star formation rate over the last lO^yr as 
determined from the stellar continuum has proven to be successful 
in reproducing star formation rates as estimated from emission 
lines ( |Asari et al.|2007 >. 



derived quantities we set the maximum of agcmass and min- 
imum of Av from |Cid Fernandes et~aL] p005) to 13.5Gyr 
and O.Omag, respectively, as opposed to 18Gyr and -0.9mag 
in the catalogue. Note, however, that agCmass derived from 
spectroscopy may still be systematically older because the 
star formation histories from which agCmass are derived in- 
clude SSPs with ages up to 20Gyr. Additionally, we set the 
mmimum SSFR to -3 in loglO(M,^n/(10^°Gyr)) as non star 
forming galaxies in our photometric model have SSFR—-ln{ 
/(lO^^Gyr)). Fig. 12 shows a comparison of 



in loglO(M,,n/(10'^Gyr)). Fig, 

derived parameters in the {i-z) vs. [g-r) diagram. The overall 
agreement between the spectroscopically and the photomet- 
rically determined quantities appears to be mediocre. The 
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Figure 11. Illustration of galaxies with special characteristics 
according to our BC03hr model with dust extinction. From all 
galaxies (black dots) objects with high expected extinction (green 
dots), low expected extinction (grey dots) and galaxies falling off 
the model grid (yellow dots) are selected. SDSS color composites 
are shown for randomly chosen examples from each group. Note 
the high frequency of highly inclined disks in the high extinction 
sample. 



main differences in the behavior in {i-z) vs. [g-r) are that: 
1) photometric metallicities increase with (i-z) and are al- 
most independent of (^-r), while spectroscopic metallicities 
also depend on the latter; 2) very young objects, as deter- 
mined from photometry, are present at (^-r,i-2;) = (1.0,0.5), 
but these are not seen if spectroscopic ages are used. A di- 
rect comparison between the derived parameters, see Fig. 
|13| does show a significant scatter as well as some system- 
atic differences. As the differences are rather "irregular" we 
refrain from an attempt to parameterize them. 

Under the assumption that the spectroscopically de- 
rived quantities are of higher accuracy than those derived 
from photometry we can probe the quality of the parame- 
ters included in the fitting. Results of such a test are shown 
in Table [2] in terms of the standard deviation in the spectro- 
scopically derived quantity at fixed photometrically derived 
values. The use of i/^,^,r,i,z-band photometry alone leads to 
a low standard deviation in all properties except the SSFR. 
This problem is solved if a prior as a function of Mr is intro- 
duced and additionally a slight improvement can be achieved 
by introducing a/6 in the fitting. Note however that Mr can 
be estimated using i^,^,r,i,z-band photometry with an ac- 
curacy of about Imag, a sufficiently good estimate to reach 
almost the same accuracy as when Mr itself is used. Us- 
ing only Mr, or Mr and a/6, we can still constrain sev- 
eral quantities as the model galaxy properties have a strong 
mass/ luminosity dependence as well as some structural de- 
pendences ( Guo et ar||2011 ), but the standard deviation is 
considerably higher in many of the parameters. 

Moreover, with a similar technique we can test whether 
the exclusion of certain bands in the fitting can improve the 
quality of the derived parameters. We find that this is the 



case. Removal of one or two bands can indeed decrease the 
standard deviation of the comparison with the spectroscopy 
by 5-10% for each of the parameters tested. However, an in- 
crease in the accuracy of one parameter typically comes at 
the expense of a diminished accuracy in another. For exam- 
ple, excluding u and g leads to a better estimate in the mass 
weighted metallicity but makes the accuracy of the SSFR 
much worse. As no overall improvement is achieved by ex- 
cluding one or two bands, we prefer to keep all bands in the 
fitting. 



6.4 Insights from spectral fitting 

While the main focus of this paper is testing model predic- 
tions for i^,^,r,z,z-band photometry an efficient test can be 
achieved by comparing observed and modelled spectra in- 
stead of colors due to the much larger information content. 
A drawback is that we already expect there will be some 
issues with the modelling because the chemical content of 
a galaxy is described by a single parameter in our model, 
i.e., the metallicity, an approximation that is not expected 
to reproduce all possible line ratios in galaxies. 

We randomly select 1000 galaxies in our sample for 
which we obtain fully reduced flux calibrated SDSS spec- 
tra. These have a median signal to noise ratio per pixel of 
19 in the r-band and are fitted to our model library in the 
wavelength range 3700A to 7000A (rest frame), as described 
in the Appendix |B] Sky lines and emission lines are masked 
out in the fit and the line of sight velocity distribution is 
taken into account. 

The median reduced of the fit to the 1000 galaxies 
is 1.4 for the BC03hr models with dust. A visual inspection 
of observations and best-fitting models shows that not only 
issues with lines, but also with the continuum fit exist. We 
note that the median reduced can be substantially low- 
ered if the fitted region is restricted to a shorter wavelength 
interval like 5000A to 5500A instead of 3700A to 7000A 



7 DISCUSSION 

Before starting the discussion about the model performance 
and prediction it is worth taking a look at how the models 
are constructed. 



7.1 Model ingredients 



As opposed to [Conroy et al.| ( |2009| '2010) and 'C onroy 



lorunn {2010} our aim has not been to judge each model in- 
gredient and we thus find it sufficient to say that some, but 
certainly not all, SSP models can be successfully used to pre- 
dict i^,^,r,i,z-band photometry of the local galaxy popula- 
tion. Caution should be taken in the choice of model. Of the 
models we tested (see Table [T]), we recommend the Bruzual 
& Chariot \200^ high resolution models for the purpose of 
predicting optical broad hand photometry of galaxies. 

Does our library of star formation histories make sense? 



^ Other spectral fitting methods that do not rely on a set of 
predefined models (in our case 50000) have been able to produce 



better fits to SDSS spectra (c.f. Cid Fernandes et al.|2005 l 
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Figure 12. Comparison between stellar mass weighted ages, age^^^^, r-band luminosity weighted ages, age^, stellar mass weighted 
metallicities, Zmass, luminosity weighted metallicities, zi^^m specific star formation rates, SSFR, r-band mass to light ratios, M/Lr, and 
total y-band extinctions. Ay, for the spectroscopic comparison sample (greyscale histogram) as derived from u,g,r,i,z-hdind photometry 
(phot) and spectra (spec) ( |Cid Fernandes et al.|2005( [Gallazzi et al.|200 5). The red lines show the 1:1 relations, blue lines indicate the 
mean in the spectroscopic values in bins of the photometric values. The bin widths are 1/30 of the entire range displayed, however, these 
are in some cases adjusted to enclose at least 50 galaxies in each bin. On top of each figure the standard deviation in the spectroscopic 
values around the mean, i.e. the blue line, is indicated. 



Table 2. Dependence of the standard deviation in stellar mass weighted age, age-mass , luminosity weighted age, agcr, stellar mass 
weighted metallicity, Zmass, luminosity weighted metallicity, ziurm specific star formation rate, SSFR, r-band mass to light ratio, M/Lr, 
and y-band extinction. Ay, compared with the spectroscopically determined values. The parameters used in the computation of the 
maximum likelihoods are w,gf,r,i,2;-band photometry, ugriz, a prior based on Mr, prior(Mr), a prior based on a Monte Carlo estimate of 
Mr from u,g,r,i,z-hand photometry, prior (ugriz— )► Mr), and the isohotal axis ratio a/b. 



rms 

fitting parameters agcmass agcr Zmass zium SSFR M/Lr Ay 

(Gyr) (Gyr) {zq) {zq) (dex) {Mq/Lq) (mag) 



ugriz 


2.8 


1.6 


0.33 


0.43 


0.87 


0.49 


0.21 


prior (Mr) 


2.9 


2.3 


0.31 


0.44 


0.50 


0.89 


0.27 


a/b + prior (Mr) 


3.2 


2.3 


0.38 


0.48 


0.73 


1.01 


0.24 


ugriz + prior (ugriz— ^ Mr) 


2.8 


1.5 


0.34 


0.44 


0.75 


0.50 


0.21 


ugriz + a/b + prior (ugriz— )► Mr) 


2.8 


1.5 


0.34 


0.44 


0.49 


0.50 


0.21 


ugriz + prior (Mr) 


2.9 


1.6 


0.32 


0.44 


0.49 


0.50 


0.21 


ugriz + a/b + prior(Mr) 


2.8 


1.5 


0.32 


0.44 


0.49 


0.50 


0.21 



Given the distribution of average values for our models it 
seems we have covered essentially all possibilities with ages 
from 0-13Gyr, metallicities from 0.00-0.05 in Z and with 
extinctions from 0-3mag. Furthermore we note that Eq. [l] 
produces star formation histories that resemble results from 
modelling of spectra of nearby galaxies by | Heavens et al.| 



(2004). As mentioned in the introduction the inclusion of 



some stochasticity is expected to improve the realism of the 
models. In the context of the amplitude of the fluctuations 
in star formation over time we note the following. By study- 
ing the ensemble of UV to Ha flux ratios in star forming 



dwarf galaxies Lee et al. (2009) investigated the starburst 
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Figure 13. Galaxy properties for the spectroscopic comparison sample as we derived from w,5f,r,i,2;-band photometry (phot) and |Cid| 
[Fernandes et al.| ( |20Q5| >; [Gallazzi et al.| ( [2005j derived from spectra (spec). Color coded pixels show the average values of stellar mass 
weighted ages, age^^^^, r-band luminosity weighted ages, age^, stellar mass weighted metallicities, 7,mass, r-band luminosity weighted 
metallicities, z^^^, specific star formation rates, SSFR, r-band mass to light ratios, M/Lr, and total V-band extinctions. Ay, in the 



(i-z) vs. (g-r) diagram. Spectroscopic quantities are derived either using spectral synthesis techniques (age, 
absorption line index analysis (zium)- 



raass-i ^rnass- 



age^) or an 



mode of dwarf galaxies. They conclude that starbursts with 
amplitudes of ^ 4 above the quiescent mode, which last for 
about lOOMyr and occur every l-2Gyr, could explain the 
whole population. The majority of our models do not have 



starburst episodes that exceed 4 times the local time average 
(see Fig.fTO). 



Our dust model is rather simple. However, given the 
success in reproducing the optical colors of the observed 
galaxy population (after some fine-tuning of constants, see 
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Fig. [6|, along with the observations that partially support 
our dust model (see Section [s]) we believe it captures some 
of the essentials of obscuration by dust. It also turns out 
that the distribution of z-band extinctions for our models 
looks rather similar to the values for the galaxy sample of 
Kauffmann et al.| (12003). 



In terms of line emission Fig. |3] suggests that modelling 
Ha and H/3 is not sufficient for the modelling of optical col- 
ors. If GALEV treats line emission adequately, this model 
further tells us that line emission only has a second order 
effect on integrated i^,^,r,i,2;-band photometry of nearby 
galaxies (Mr ^ 17), but caution should be taken at higher 
redshifts since migration of strong emission lines between fil- 
ter bands does cause color shifts. Future studies improving 
the accuracy of the estimation of emission line strengths are 
on the way ( [Gyory et al.|201ip and will further improve the 
quality of stellar population models. 

Given the limited sensitivity of i/^,^,r,i,z-band photom- 
etry to the IMF as shown in Fig. [2] we conclude that 
the IMF slope, and consequently the stellar mass, cannot 
be measured with integrated optical broad band colors of 
galaxies at least in the range explored here (single slope 
IMFs from 2.1-2.9, Salpeter=2.35). This result is expected 
considering that most light from a stellar population with 
a "normal" IMF comes from stars of a rather limited mass 
range that reside close to the main sequence turnoff, on the 
giant branches or that are supergiants. 

Is the number of models sufficiently large? When de- 
termining galaxy properties along with errors from colors 
using, e.g., a Bayesian maximum likelihood, the density of 
models in color space is of importance. Judging from Fig. 
|7|9| our library is sufficiently large for this purpose except at 
the very edges of the grid if using (u-g) vs. (r-i) or (g-r) vs. 
(i-z). However, if more than two colors are used the required 
density of the model library increases and caution has to be 
taken. The errors in Fig. ^] do not appear to be underesti- 
mated, except for one galaxy falling outside the model grid, 
and the method we use thus seems to work properly even 
for 10 colors. 

The galaxy spectra are rather poorly reproduced with 
our model library, as mentioned in Section [6^ both due to is- 
sues with lines and the continuum. The fact that the reduced 
can be lowered by considering a shorter wavelength range 
probably indicates that our model struggles to reproduce 
the overall continuum shape, but can preform better over 
shorter wavelength intervals. Spectral lines may be poorly 
reproduced due to the simple one parameter treatment of 
element abundances. We should on the other hand be able 
to adequately model the continuum, especially considering 
its close connection to colors. However, in the modelling of 
galaxy spectra it is often assumed that the continuum is un- 
reliable due to, e.g., flux calibration issues, and continuum 
variations are therefore removed with polynomial fits (cf. 
Chihngarian et al.|[2007| |Koleva et"aLl|2009j ) . It is thus not 
completely clear whether the poorly fitted spectra mainly 
refiect problems with the models or the observations. 



7.2 Model predictions 

It is well known that optical colors can be used to separate 
star forming galaxies from quiescent ones. However, metal- 
licity and dust extinction also have major impact. Our mod- 



els suggest that, by using optical colors only, it is possible 
to reduce the age-metallicity degeneracy and constrain ages 
and metallicities fairly well over a surprisingly large part of 
the parameter space (see Fig. 7][8 ). Furthermore, an inter- 
esting prediction of the model is that, to some extent, the 
amount of internal extinction can be measured using optical 
colors (Fig. [9]). Is this really the case? We expect extinc- 
tion to have the strongest effect in edge on disk galaxies. 
We therefore selected two subsamples of galaxies, with high 



and low estimated extinction, respectively (Fig. 11). It does 



indeed look like the high extinction bin contains a lot more 
edge on disk galaxies than the low extinction bin and, more- 
over, dust is seen in many more cases in the high extinction 
bin. It thus appears as if our color based estimate is reason- 
able. 

Fig. |12| shows a rather good agreement between spectro- 
scopically and photometrically determined galaxy properties 
in terms of average ages, metallicities, effective extinctions, 
specific star formation rates, and mass to light ratios. The 
strong correlation between (i-z) and metallicity in the pho- 
tometric model (Fig, [s] and is not as pronounced in (i-z) 
vs. spectroscopic metallicity (Fig. 12), suggesting that this 



particular model prediction may not be robust. We note, 
however, that the correlation between (i-z) and metallicity 
is present in all models we have considered (see Table[T]). The 
wavelength dependence of the luminosity weighting, i.e. the 
wavelength dependence of the relative contribution of light 
from young and old stellar populations, could infiuence the 
metallicity estimate. If the spectroscopic metallicities are 
based on strong absorption lines line like e.g. the Mg line 
at ^ 5IOOA, they are likely more infiuenced by young stel- 
lar populations and blue horizontal branch stars than the 
i and ^-band which have effective wavelengths of 7472 and 
8917A. The source of the discrepancy in metallicity is thus 
not resolved. 

As the amount of dust present is coupled to a galaxy's 
specific star formation rate in our model, galaxies with 
young stellar populations are often heavily obscured. The 



introduction of the prior (see Sect. 6.1 ) reduces the degener- 
acy between old and dust-free and young and dusty galaxies. 
However, a small number of old galaxies (ages determined by 
spectroscopy) still have young stellar populations according 
to our photometric model (Fig. |12| ). Figure 13 reveals some 
systematic differences between the parameters derived from 
photometry and spectroscopy along with some scatter. How- 
ever, deviations are to be expected given the large difference 
in the amount of information contained in the two kinds of 
data (5 data points along the spectral energy distribution for 
u,g,r,i,z with respect to ^2000 for SDSS spectra). To probe 
galaxy properties with optical photometry thus stands as a 
viable option, especially useful when limited spectral infor- 
mation is available such as in the case of severe aperture 
bias. We would also like to point out that spectroscopically 
derived quantities are not necessarily "correct" as model de- 
pendencies may be present. 

The exclusion of some bands can slightly diminish the 
differences between photometrically and spectroscopically 
derived quantities (see Sec. 6.4). This may be caused by 



the wavelength dependent luminosity weighting, but could 
also refiect that certain colors do not trace all parameters 
we study. 

An interesting result is that the newer CB07 models per- 
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form worse than the BC03 models in reproducing the colors 
of the nearby galaxy population. We would like to point out 
that this is not the same as saying that the ingredients of 
the BC03 model are somehow better than the newer version. 
We speculate that the offsets in the CB07 models may be 
a backlash of trying to optimize the model for predictions 
over a larger wavelength range. In BC03 as well as CB07, 
models based on the empirical library of stellar spectra by 



Le Borgne et al. (2003) do improve the modelling with re- 



spect to models based on stellar atmosphere models of |Leje-| 
une et al. (1997). This suggest that more uncertainties are 
introduced through a complete modelling than through an 
empirical approach (see also Maraston et al.|2009| [Maraston 



StrQmback||2011| [Peacock et al.||2011| ) 

A subclass of galaxies at (^-r,z-z) = (0.4,0.0) can not 
be modelled by any of the models considered (see Fig. 
11). Most of these galaxies are rather faint in the z-band 
(< niz 17.0 as compared to < niz >~ 15.5 for the entire 
sample) , which has led us to suspect that a part of this popu- 
lation could simply be an artifact caused by photometric er- 
rors. The fraction of galaxies with small z-band photometric 
errors (Jm, < 0.05 in the region (0.3 < (g-r) < 0.4,-0.05 < 
(i — z) < 0.05) is a factor of two smaller than for the overall 
sample. However, this is not sufficient to explain the entire 
population. Another explanation is that these galaxies are 
even more metal poor than the most metal poor models 
(see Fig.jsj) or, alternatively, that the SSP models fail at low 
metallicity. 

The star formation histories, along with average ages, 
metallicities and r-band optical depth shown in Fig. [lO] 
look reasonable at first glance. A quiescently looking galaxy, 
which in fact does not have any emission lines in its SDSS 
spectra, should according to the models have some amount 
of recent star formation. However, non star forming models 
lie within the errors. We thus conclude that the model out- 
put is reasonably reliable - within the limits of broad band 
photometry - in providing information about a galaxy's star 
formation history and dust content. 



8 SUMMARY 

If an adequate range of star formation histories, chemical en- 
richments and dust obscuration are being considered, some 
of the currently available stellar population models perform 
well in reproducing the integrated optical u,g,r,i,z-hsind 
photometry of nearby galaxies (see Fig.[6}. However, strong 
differences between the various models are seen and it is 
therefore necessary to carefully consider the choice of model 
before proceeding with the analysis. Our preferred model li- 
brary is available on request. Optical broad-band colors are 
insensitive to the slope of the IMF, but can on the other 
hand put constraints on star formation histories, chemical 
enrichments, and dust extinctions as verified through a com- 
parison with spectroscopic data. The accuracy of the con- 
straint depends on the galaxy's position in color space. This 
is mainly due to overlap of reddened and dust-free models 
in certain parts of the color space. 
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APPENDIX A: GALAXY COLORS 

We test the colors of the galaxy sample we use in this paper 
against two independent samples derived from SDSS data. 



The first sample from Lisker et al. (2008); Janz & Lisker 
([2009) and Meyer et al. (in prep.), hereafter SI, consists of 
Virgo cluster galaxies for which colors are measured within 
2 effective radii or 1 Petrosian radius. The second sample 
is from Hansson et al. (in prep)., hereafter S2, and consists 
of an essentially volume limited sample of galaxies brighter 
than Mr — 16 within 50Mpc. Colors for these galaxies are 
measured within 2.5Kron radii, as determined from the r- 
band image using SourceExtractor CBertin & Arnouts 1996 ). 



For these two samples as well as for the sample we use in this 
paper, hereafter S3, (r-i) vs. (u-g) and (i-z) vs. (g-r) colors 
are given in Fig. |A1| The general agreement is good, but 
differences are seen. S3 have on average higher values of {i- 
z) at [g-r)— 0.7. We test whether this can be explained by the 
difference in how the magnitudes are computed by looking 
at the change in colors when adopting Petrosian magnitudes 
(PetroMag) instead of modelMags. The median change in {i- 
z) is -0.06mag and for the other colors it is between -0.02mag 
and O.OOmag. Because the definition of Kron magnitudes is 
similar to Petrosian magnitudes we conclude that this offset 
can be explained by the difference in magnitude definitions. 
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37OOA to 7OOOA. Regions affected by the four brightest sky- 
Unes (015577, OI6300, 016354 and NaD5890) are masked 
out. The model hbrary is at a sUghtly lower resolution 
(cFmod = 3A/2.35) than the data (aobs ^ 2.5A/2.35 at 
A = 5OOOA) and we therefore degrade the resolution of the 
observed spectra to FWHM=3A by convolving with a Gaus- 
sian with a — \J cr'^^^ — crl^^. As the model library is fairly 
large a direct fit to all models would be rather computa- 
tionally expensive. We therefore carry out a two step fitting 
making use of binned models produced in the following way. 
The rms between one model and all other models is com- 
puted. The average of the 225 model with the lowest rms 
is then computed, constituting the first binned model. The 
procedure is repeated but all models already included in the 
binned spectra are excluded. After 223 repetitions we have 
223 binned models, each of which, except the last one, points 
to 225 other models. In the fitting procedure we first find the 
best binned model and secondly the best model in the bin. 
Such a fitting need not be equivalent to a of ^lH mod- 

els since the best-fitting model not necessarily needs to be 
contained in the best-fitting bin. Nevertheless, we find this 
approach preferable since it accelerates the fitting procedure 
by a factor of about 10^. 

Two things need to be considered before comparing the 
models with the observations, namely the line of sight ve- 
locity distributions of the stars and the emission lines. We 
assume that the former can be modelled by a Gaussian dis- 
tribution. For each spectrum the minimum is found after 
convolving the spectrum with Gaussians having dispersions 
from 100 to 285 km/j^in 12 logarithmically spaced steps. 
The procedure is repeated 25 times after having perturbed 
each pixel of the observed spectra with its errors. The best- 
fitting dispersion is taken to be 



Figure Al. SDSS pipeline modelMags for the galaxy sample 
used in this paper (grey 2D histogram), Kron magnitudes of SDSS 
galaxies from Hansson et al. (in prep.) and Petrosian magnitudes 
from ( [Lisker et al.||2008| |Janz Lisker||2009| Meyer et al. (in 
prep.)) in the (r-i) versus (u-g) and (i-z) versus (g-r) planes. All 
magnitudes were corrected for Galactic foreground extinction. 



SI also have many more galaxies in the region around (g- 
r,i-z) = (0.6,0.0). We find that all these galaxies are fainter 
than about Mr = — 16 and as such are therefore absent in 
S3. 



Y^aj exp(-|x?) 
Eexp(-|x?) 



(B2) 



where ai is the best-fitting dispersion of the ith Monte 
Carlo realization, and where Xi is the lowest reduced of 
the ith Monte Carlo realization. 

The emission lines we consider in the fit are those 20 
emission lines considered in Stoughto n et al.| (12002 ) . These 
lines are expected to be the strongest in star forming galaxies 
according to the models of [Anders Fritze-v. Alvensleben] 
( (2003| ) , and furthermore include the strongest emission lines 
expected due to AGN activity. For simplicity all 20 lines are 
masked out in the fit. 



APPENDIX B: SPECTRAL FITTING 

The spectra are deredshifted using the redshift provided by 
the SDSS according to 



■ v/c 



1 +v/c 



(Bl) 



where Xrest and Aobs are the rest wavelength and observed 
wavelength, respectively, v is the velocity corresponding to 
the Doppler shift and c is the speed of light. 

The best models in our library of BC03hr spectra re- 
producing the data are found by means of calculating the 
reduced between each model and the observations. In 
the fitting process we make use of the spectral region from 



APPENDIX C: BAYESIAN MAXIMUM 
LIKELIHOOD MODELLING 

Assuming the errors are Gaussian the star formation his- 
tory, stellar mass weighted stellar metallicity and effective 



The number of models in each bin was chosen for speed opti- 
mization (VSOOOO ^ 224). 

-'^-'^ A lower limit of lOOkm/s is adopted to avoid ending up with 
observations at higher resolution than the models at short wave- 
lengths (SDSS spectra have 6X/X ^ 2000). With an upper limit 
of 285km/s we can model all except some of the most massive 
galaxies in the Universe ( Bernardi et al.|2007| >. 
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extinction can be computed according to 



-24.5<M <=-22.5 



-23.5<M <=-21.5 



i—nj = m — -(■ 

i=l j = l 



i=n j—m 
1=1 3 = 1 



(Cl) 



where a is any of the properties mentioned above, i denotes 
a model, o an observation, Cj a color, (Tcoj the error in an ob- 
served color and the summations are done over all n models 
and all m colors. The prior, i.e. the weights, Wi^ are taken as 
the ratio of the numbe r of models in the library and the num- 
ber of galaxies in the |Guo et al. (2011) catalo for binj^ 
in Mr where the ratios are computed for models and galax- 
ies within 0.1 dex in stellar mass weighted age and 0.3dex in 
stellar mass weighted z of i. Figure |C1| illustrates how the 
Bayesian approach influence the model grid in terms of age 
and metallicity. 
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